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SUMMARY 

The  orbit  of  Cosmos  268  rocket  (1969-20B)  has  been  determined  at  28  epochs 
during  its  342-day  life,  with  the  aid  of  the  PR0P5  orbit  refinement  program. 

All  available  observations  were  used,  including  16  from  the  Hewitt  camera  at 
Malvern,  28  from  the  200-mm  camera  at  Meudon,  56  from  the  kinetheodolite  at  the 
Cape  Observatory,  700  visual  observations  from  volunteer  observers,  500  US  Navy 
observations  and  200  British  radar  observations.  The  orbits  are  of  very  good 
accuracy  for  such  a  high-drag  satellite,  most  of  the  values  of  inclination 
having  standard  deviations  less  than  0.002°.  The  most  accurate  orbits  are 
those  utilizing  photographic  observations,  and  the  best  of  these  has  standard 
deviations  of  0.00001  in  eccentricity  and  0.0001°  in  inclination. 

The  values  of  inclination  obtained,  after  correction  to  allow  for  the 
effects  of  other  perturbing  forces,  have  been  analysed  to  determine  zonal  wind 
speeds  in  the  upper  atmosphere  at.  heights  a  little  above  perigee  (230-250  km) 
averaged  over  latitudes  up  to  about  25°.  The  results  show  a  clear  distinction 
between  the  wind  at  night  (21  h  to  03  h  local  time) ,  which  is  west-to-east 
with  an  average  speed  of  140  t  50  ra/s,  and  the  wind  by  day  (04  h  to  19  h), 
which  is  east-to-west  with  an  average  speed  of  110  ±  50  m/s. 


Departmental  Reference:  Space  381 


2 


202 


CONTENTS 

Page 

1  INTRODUCTION  3 

2  THE  OBSERVATIONS  3 

3  THE  ORBITS  4 

3.1  General  review  4 

3.2  Techniques  for  orbit  improvement  5 

3.3  Perigee  height  6 

3.4  Decay  rate  7 

4  UPPER-ATMOSPHERE  WIND  SPEEDS,  FROM  THE  CHANGES  IN  INCLINATION  7 

5  RESIDUALS  OF  THE  OBSERVATIONS  10 

6  CONCLUSIONS  12 

Table  1  Orbital  parameters  for  Cosmos  268  rocket  (1969-20B) 

with  standard  deviations  13 

Table  2  Values  of  A  obtained  9 

Table  3  Analysis  of  residuals  for  selected  stations  11 

References  15 

Illustrations  Figures  1-4 

Detachable  abstract  cards 


Drgs:  C04/903312  to  903315 


202 


3 


1  INTRODUCTION 

Cosmos  268  was  launched  on  5  March  1969  into  an  orbit  inclined  at  48.4° 
to  the  equator,  with  an  initial  orbital  period  of  109.1  minutes.  Its  final- 
stage  rocket,  designated  1969-20B,  entered  a  similar  orbit,  with  an  initial 
perigee  height  220  km  and  apogee  height  nearly  2200  km.  This  rocket,  8  m  long 
and  1.65  m  in  diameter,  was  a  good  object  for  optical  observation  from  Britain, 
its  stellar  magnitude  at  45°  elevation  being  about  + 2  near  perigee  and  +5  near 
apogee.  The  orbit  contracted  rapidly  because  perigee  was  so  low,  and 
1969-20B  decayed  on  11  February  1970.  During  its  342-day  life  the  orbital 
period  decreased  by  more  than  20  minutes,  so  the  decrease  in  orbital  inclina¬ 
tion  caused  by  atmospheric  rotation  was  expected  to  be  quite  substantial, 
perhaps  as  much  as  0.1°.  Consequently  the  orbit  seemed  promising  for  use  in 
determining  the  upper-atmosphere  rotation  rate,  and  1969-20B  was  given  a  high 
priority  for  observing  by  the  British  optical  and  radar  tracking  stations, 
including  the  Hewitt  camera  at  Malvern. 

The  orbit  of  1969-20B  has  been  determined  from  all  available  observations 
with  the  aid  of  the  RAE  orbit  determination  program  PROP*  in  the  PR0P5  v  ;rsion, 
at  28  epochs  between  5  April  1969  and  10  February  1970.  This  Report  describes 
the  orbit  determinations,  and  analyses  the  resulting  values  of  inclination  to 
determine  atmospheric  rotation  rates  at  heights  near  240  km. 

2  THE  OBSERVATIONS 

Six  groups  of  observations  were  available.  The  most  accurate  were  those 
from  the  600-mm  Hewitt  camera  at  Malvern.  Although  these  observations  do  not 
exploit  the  full  accuracy  of  the  camera,  because  the  plates  are  reduced  by  a 
quicker  method  than  is  used  for  geodetic  work,  the  Hewitt  camera  observations 
are  nevertheless  usually  accurate  to  about  2  seconds  of  arc.  For  1969-20B 
there  were  16  of  these  observations,  from  4  plates. 

The  second  group,  almost  as  accurate,  were  the  photographic  observations 
from  the  200-mm  camera  at  Meudon.  There  were  28  of  these  observations,  on 
7  transits,  with  an  accuracy  of  about  5  seconds  of  arc. 

The  next  most  accurate  group  of  observations  was  that  from  the  kine- 
theodolite  at  the  Cape  Observatory,  South  Africa,  with  accuracy  of  about 
li  minutes  of  arc.  These  observations,  56  in  number,  were  of  particular 
importance  because  of  their  geographical  location:  since  northern-hemisphere 
observations  are  usually  predominant,  the  addition  of  a  few  accurate  southern- 
hemisphere  observations  can  greatly  improve  the  orbital  accuracy. 
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Observations  from  the  other  three  sources,  which  made  up  the  great 
majority  of  the  data,  were  all  of  rather  similar  accuracy,  between  lj  and 
5  minutes  of  arc  in  direction  and  (for  radar  only)  between  1  and  2  km  in 
range.  The  most  numerous  of  these  (over  700)  were  visual  observations  from 
the  volunteer  observers  who  report  to  the  Radio  and  Space  Research  Station, 
Slough  or  the  Moonwatch  Division  of  the  Smithsonian  Astrophysical  Observatory, 
and  from  the  theodolites  at  Meudon  and  associated  stations  in  France.  Since 
most  of  the  visual  observations  were  made  from  latitudes  between  48°  and  52°, 
they  had  a  strong  influence  in  fixing  the  value  of  the  orbital  inclination, 
on  those  occasions  when  camera  observations  were  not  available. 

Another  major  source  of  observations  was  US  Navy  data,  kindly  supplied 
by  the  US  Naval  Research  Laboratory.  More  than  500  US  Navy  observations  were 
used,  and  many  of  the  orbit  determinations  would  have  been  impossible  without 
them. 

The  sixth  category  of  observations  comprised  those  from  British  radar 
stations,  particularly  the  radar  tracker  at  RRE,  Malvern.  These  observations 
totalled  nearly  200  and  were  particularly  useful  at  the  end  of  the  satellite's 
life. 

3  THE  ORBJVS 
3.1  General  review 

The  28  orbits  obtained  are  given  in  Table  1  on  page  13.  As  usual  with 
PROP,  the  epoch  is  at  0  h  on  the  chosen  day  and  the  definitions  of  the  main 
orbital  elements  are  given  below  the  table.  In  PROP  the  mean  anomaly  M 
is  fitted  by  a  polynomial  of  the  form 

M  =  Mq  +  Mxt  +  M2t2  +  M3t3  +  M4t4  +  M5t5  (1) 

where  t  is  the  time,  '.assured  from  the  chosen  epo'h,  and  the  number  of 
coefficients  M2>  M^,  ...  used  is  open  to  choice.  Wi,,h  a  high-drag  orbit 
like  that  of  1969-20B,  the  best  number  of  coefficients  to  use  has  to  be 
found  by  trial  and  error.  For  4  of  the  runs  the  best  results  were  obtained 
by  using  M^  to  M2  only;  13  of  the  runs  required  MQ  to  M^;  10  required 
Mq  to  M^;  and  for  the  last  run,  near  the  end  of  the  life,  the  full  comple¬ 
ment  of  Mq  to  Mj  was  found  to  be  necessary.  The  orbits  given  in  Table  1 
fit  the  observation.'  well,  with  e,  the  parameter  measuring  the  goodness  of 
fit,  usually  rather  less  than  1. 
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The  7  orbits  which  utilised  Hewitt  or  Meudon  camera  observations  are 

marked  with  asterisks  in  Table  1,  and  these  can  appropriately  be  discussed 

first.  For  one  of  these  orbits.  No. 5,  the  supporting  observations  were 

unfortunately  quite  inadequate  both  in  number  and  in  geographical  coverage, 

and  the  orbit  obtained  was  rather  mediocre.  The  remaining  6  orbits  gave  values 

of  inclination  which  were  of  excellent  accuracy,  with  the  standard  deviations 

ranging  from  0.0001°  and  0.0007°,  corresponding  to  distances  between  10  and 

90  metres.  One  orbit.  No. 14,  has  the  lowest  standard  deviations  yet  obtained 

2 

for  such  a  high-drag  orbit  (Mj  =1.3  deg/day  ),  with  sd  in  eccentricity  of 
0.00001,  corresponding  to  sd  in  perigee  height  of  70  m,  and  sd  in  inclination 
of  0.0001°.  Orbit  No. 7,  which  utilizes  4  plates  from  the  Meudon  camera  all 
within  a  time  interval  of  48  hours,  has  even  lower  sd  for  some  of  the  para¬ 
meters,  0.00006°  for  i,  0.001°  for  0  and  0.0002  deg/day  for  M^;  but  the 
observations  are  few  in  number  (25)  and  are  all  from  latitudes  between  48° 
and  54°,  so  that  there  is  more  danger  of  bias  which  may  considerably  exceed 
the  sd.  Also,  as  would  be  expected  with  such  limited  orbital  coverage,  the  sd 
in  eccentricity  is  relatively  large  (0.00006). 

The  21  orbits  obtained  without  the  aid  of  camera  observations  were  of 
somewhat  varied  accuracy,  depending  on  the  number  and  geographical  location 
of  the  observations.  On  14  of  these  orbits  (including  the  final  one)  the 
sd  in  inclination  was  between  0.0013°  and  0.0020°,  while  the  other  7  had 
values  between  0.0022°  and  0.0038°.  Standard  deviations  in  eccentricity 
range  between  0.00001  and  0.00006  on  all  the  orbits  except  run  9,  for  which 
the  geographical  coverage  was  very  poor.  The  sd  in  the  right  ascension  of 
the  node  was  between  0.002°  and  0.007°,  again  with  the  exception  of  run  9. 

3.2  Techniques  for  orbit  improvement 

It  was  usually  necessary  to  run  the  PROP  program  several  times  at  each 
epoch,  wit  .  each  run  involving  a  number  of  iterations  in  the  refinement 
process,  before  arriving  at  the  final  orbits  listed  in  Table  1.  Quite  often 
the  first  attempt  would  end  in  divergence,  because  the  initial  orbit,  usually 
predicted  from  the  previous  epoch,  was  not  accurate  enough.  On  average,  the 
first  successful  attempt  needed  about  10  iterations,  the  final  convergence 
often  being  delayed  because  of  the  presence  of  observations  which  were  in 
error  by  a  fairly  small  amount  and  tended  to  be  rejected  one  by  one  in  succes¬ 
sive  iterations.  Usually  the  first  attempts  were  made  with  the  parameters  MQ 
to  in  the  polynomial  for  M,  but  it  was  often  necessary  to  re-run  with 
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Mg  to  M2  only,  or  with  M^  to  M^,  to  see  whether  the  fit  was  better  with 
these  sets  of  coefficients.  The  other  main  reason  for  re-runs  was  the  survival 
of  observations  which  'looked  wrong',  but  were  not  bad  enough  to  be  rejected 
by  the  computer .  Observations  are  rejected  if  they  have  a  residual  greater 
than  3e  times  the  assumed  error:  but  if,  say,  5  out  of  30  observations 
have  weighted  residual  3,  their  initial  acceptance  inevitably  keeps  e  greater 
than  /45/30,  i.e.  greater  than  1.2,  thereby  ensuring  their  survival. 

When  the  'final'  orbits  had  been  arrived  at,  and  the  values  of  inclination 
were  plotted  against  time,  it  appeared  that  some  values  near  the  end  of  the  life 
were  out  of  step  with  the  general  trend.  These  orbits  were  re-examined  and 
re-run  with  various  modifications:  the  orbit  eventually  fitted  the  observations 
better  -  and  better  values  of  inclination  were  obtained  -  usually  after  omit¬ 
ting  observations,  either  at  the  beginning  or  end  of  the  set  used,  or  at  a  time 
when  irregularities  in  drag  seemed  to  be  spoiling  the  fit.  All  these  diffi¬ 
culties  are  most  probably  due  to  the  polynomial  for  M  failing  to  fit  the  real 
drag,  which  was  not  only  irregular  but  also  very  large,  with  M2  greater  than 
2  dag/day . 

3 . 3  Perigee  height 

Fig.l  gives  an  impression  of  the  variations  in  perigee  height..  The 
crosses  indicate  values  of  a(l  -  e)  -  R,  i.e.  the  perigee  'height'  with 
respect  to  the  Earth's  equatorial  radius  R  (=  6378.2  km).  The  broken  curve 
drawn  through  the  crosses  shows  the  oscillation  in  perigee  distance  produced 
by  the  odd  zonal  harmonics  of  the  Earth's  gravitational  field.  For  a  satellite 
with  inclination  48.4°  the  amplitude  of  this  oscillation^  is  5.2  km.  So,  if 
we  define  a  perigee  height  parameter,  Q,  by  the  equation 

Q  =  a(l  -  e)  -  R  +  5.2  sin  10  km  (2), 

Q  should  be  constant  apart  from  luni-solar  perturbations  and  the  gradual 
decrease  in  perigee  height  due  to  the  effect  of  air  drag.  In  Fig.l  the  values 
of  Q  for  the  28  orbits  are  indicated  by  circles,  with  sd  when  these  are 
large  enough  to  be  visible:  the  values  of  Q  are  satisfactorily  close  to  a 
smooth  curve,  though  the  rather  steep  initial  slope  of  the  curve  is  a  little 
unexpected. 

For  a  satellite  with  inclination  48°,  like  1969-20B,  the  average  Earth 
radius  over  a  complete  revolution  is  about  6372  km,  so  the  average  perigee 
height  is  about  Q  +  f  km. 
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3.4  Decay  rate 

The  decay  rate  of  1 969-20B  due  to  air  drag  was  rapid  and  variable,  and 
Tig. 2  gives  values  of 

iM  =  M2  +  3M3t  +  6M4t2  (3), 


where  t  is  the  time  from  epoch,  over  the  time  intervals  covered  by  the  obser¬ 
vations  used  at  each  epoch.  Fig. 2  indicates  the  variability  in  drag  from  week 
to  week,  but  does  not  show  the  day-to-day  irregularities,  which  tend  vo  be 
smoothed  out  by  fitting  a  polynomial  for  M. 

The  two  most  severe  magnetic  ;torms  during  the  lifetime  of  the  satellite 
occurred  on  15  May  and  30  September  1969,  and  are  marked  in  Fig. 2  by  double¬ 
headed  arrows.  Sharp  increases  in  air  density  occur  at  the  times  of  such 
magnetic  storms  and  corresponding  increases  would  be  expected  ir  M,  which  is 
directly  proportional  to  the  air  density  at  heights  near  perigei.  Fig. 2  con¬ 
firms  that  there  are  steep  increases  in  M  at  these  rimes:  the  two  orbits 
affected  (Nos. 6  and  14)  have  the  highest  (positive)  values  of  for  any 

orbit  up  to  No. 23. 

4  UPPER-ATMOSl HERE  WIND  SPEEDS,  FROM  THE  CHANGES  IN  INCLINATION 


Because  the  upper  atmosphere  is  rotating,  at  approximately  the  same  rat  i 
as  the  Earth,  the  aerodynamic  force  acting  on  a  satellite  has  a  component 
perpendicular  to  the  orbit,  which  has  the  effect  of  reducing  the  inclination 
of  the  orbit  to  the  equator  as  time  goes  >n.  For  a  high -drag  satellite, 
atmospheric  rotation  is  the  most  important  force  perturbing  i,  ant  if 
corrections  are  made  to  allow  for  lutii-solir,  odd  zonal  harmonic  and  tessera’ 
harmonic  perturbations,  the  changes  In  i  give  the  atmospheric  rotation  rat* 
and  hence  the  vind  speed  in  the  region  near  the  satellite's  perigee. 

Fig. 3  shows  the  values  of  i  from  the  28  orbits  in  Table  1,  with  -'d. 

These  values  were  corrected,  using  the  PROD  computer  program" ,  for  the  effects 

of  zonal  harmonic  and  luni-solar  perturbations,  which  vtfere  calculated  at  1-day 

intervals;  and  also  corrected  for  the  tesseral  harmonic  perturbations, 

o  ~yi  . 

which  were  always  less  than  0.0015  .  The  corrected  values  of  inclination, 
with  sd,  are  plotted  in  Fig. 4.  The  theoretical  curves  tor  the  change  in 
inclination  d"e  to  an  atmosphere  rotating  at  various  constant  rates 
(0.7,  0.8  ....  1.4  rev/day)  were  then  calculated  by  numerical  integration  of 
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the  equation  for  —  at  intervals  of  45°  in  u>  (approximately  10  days  in 
time)  initially,  and  at  closer  intervals  later  in  the  satellite's  life. 

It  immediately  became  apparent  that  the  values  of  i  in  Fig. 4  could  not 
be  fitted  satisfactorily  if  the  rotation  rate  A  were  kept  constant:  A  had 
to  be  greater  in  the  middle  of  the  life  than  at  the  beginning  or  end.  So  the 
theoretical  curve  was  divid’d  into  3  sections,  with  breaks  at  what  appeared 
to  be  the  most  suitable  dates,  namely  MJD  40431  and  40500.  At  the  first  of 
these  dates  perigee  was  at  southern  apex  (w  *  270°),  and  ~  -  0,  so  the  curve 
continues  without  any  change  in  slope.  The  second  break  occurs  when  perigee 
was  on  the  equator,  and  the  inclination  was  decreasing  most  rapidly.  The 
values  of  A  (rev/day)  which  gave  the  best  fit  for  the  combined  3-part  curve 
were  judged  to  be:  (1)  A  «  0.8;  (2)  A  *  1.3;  and  (3)  A  *  0.7.  The  estimated 

errors  in  A  are  0.1  (sd)  for  sections  (1)  and  (2),  and  0.15  for  section  (3), 
the  error  in  section  (2)  being  proportionately  smallest  -  8%,  as  against  12% 
for  section  (1)  -  because  both  its  ends  are  'fixed"  by  having  to  link  up 
with  the  other  sections. 

Since  it  is  quite  likely  that  A  varies  with  local  time,  a  scale  of 
lo«,al  c'.ae  is  given  at  the  top  of  Fig. 4.  It  is  not  possible  to  specify  an 

ex.sct  !ocal  time  for  which  A  applies,  because  (a)  the  effects  of  atmospheric 

rotation  are  appi  eciable  over  an  arc  of  about  30°  on  each  side  of  perigee 
(i.e.  ±lj  hours  approximately),  and  (b)  the  effects  are  greatest  wl en  perigee 
is  near  the  equator  and  the  curve  in  Fig. 4  is  at  one  of  its  points  of  maximum 
slope.  Thus,  to  clarify  the  picture  by  somewhat  oversimplifying  it,  the  air 
is  being  sampled  at  iui_?l  times  about  z2  hours  from  the  local  time  at  perigee 
at  the  equator  crossing;  and,  by  rather  good  luck,  the  successive  equator 
crossings  are  at  intervals  of  approximately  4  hours  in  local  time,  so  that 
the  whole  spectrum  of  local  time  is  covered,  though  with  some  bias  towards  the 
times  when  perigee  crosses  the  equator.  The  scale  at  the  top  of  Fig. 4  indicates 
that  in  section  (1)  A  applies  for  local  times  between  04  and  14  h;  in 

section  (2)  between  21  h  through  midnight  to  03  h;  and  in  section  (3)  between 

09  and  19  h*.  These  results  apply  at  a  height  J  of  a  scale  height  above  the 
mean  perigee  height  of  the  satellite:  from  Fig.l  the  mean  perigee  height 
(Q  +  6  km)  decreases  from  about  225  km  in  section  (1)  to  about  210  km  in 


*The  point  at  KJD  40503  can  conveniently  be  regarded  as  the  end  o'  section  (2), 
with  the  point  at  MJD  40522  as  the  first  in  section  (3) . 
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section  (3)  and  the  scale  height  decreases  correspondingly  from  about  32  km 
to  about  29  km.  These  results  are  summarized  in  Table  2  below. 


Table  2 

Values  of  A  obtained 


Dates 

Local  time 
at  perigee 

A 

rev/day 

Height 

km 

(1) 

1969  Apr  5  -  July  29 

04-14  h 

0.8  ±  0.1 

250 

(2) 

1969  July  29  -  Oct  6 

21-03  h 

1.3  ±0.1 

245 

(3) 

1969  Oct  6  -  1970  Feb  10 

09-19  h 

0.7  ±0.15 

230 

These  results,  which  may  be  regarded  as  averages  over  latitudes  up  to 

ii,  that  is  up  to  24°,  strongly  support  the  idea  that  the  prevailing  winds 

depend  on  local  time.  If  the  values  in  sections  (1)  and  (3)  are  combined  to 

give  A  *  0.75  ±  0.1,  there  is  a  clear  distinction  between  strong  west-to- 

east  winds  averaging  about  140  ±  50  m/s  at  night  (21  h  to  03  h)  and  east-to- 

west  winds  averaging  about  110  ±  50  m/s  by  day  (04  to  19  h).  This  result  is 

5— 8 

in  accord  with  the  somewhat  limited  evidence  from  vapour-trail  measurements 

9  10 

and  with  previous  orbital  analyses  ’  .  Most  theoretical  studies  also  suggest 

that  in  near-equatorial  regions  the  strongest  west-to-east  winds  are  likely 
to  occur  in  the  late  evening,  between  20  h  and  02  h  (see,  for  exau'nle,  Ref. 11). 

Hiller’s  orbit^of  Cosmos  307  rocket  shows  that  the  changes  ii  wind 

12 

direction  occurred  soon  after  the  two  most  severe  magnetic  storms  .  This  may 
be  coincidence,  and  it  may  also  be  by  chance  that  the  cwo  strongest  geomagnetic 
disturbances  during  the  lifetime  of  1969-20B  occurred  at  the  beginning  of 
sections  where  A  was  low.  Obviously,  the  evidence  is  as  yet  far  too  tenuous 

to  draw  conclusions,  but  the  link  between  geomagnetic  disturbances  and  zonal 

.  .  13 

wind  speeds  in  the  auroral  zone,  demonstrated  by  Rees  ,  may  lead  to  some 

indirect  influence  at  lower  latitudes. 

Since  the  local  time  at  perigee  in  Fig. 4  ranges  cvet  more  than  24  h,  it 
is  possible  to  give  a  mean  value  of  A  averaged  over  24  h  of  local  time.  The 
ob’.ious  choice  of  end  points  for  this  calculation  is  MJD  40334  and  40565:  the 
equator  crossings  covered  are  then  those  occurring  at  about  17  h,  21  h,  01  h, 

05  h  and  10  h  local  time.  Between  MJD  40334  and  40565  the  inclination 
decreases  by  the  same  amount  as  in  Fig. 4  (0.034°)  if  A  is  kept  constant  at 
a  value  of  0.9,  with  sd  estimated  as  0.1.  This  is  the  lowest  mean  value  of 
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A  so  far  obtained  from  ^bital  analysis  at  heights  between  200  and  400  km, 
and  is  significantly  lover  than  the  value  given  by  the  mean  curve  in  Ref .12. 

It  will  be  seen  that  the  last  four  points  in  Fig. 4  are  appreciably  above 

the  curve.  This  may  be  due  to  the  iifficulty  in  determining  accurate  orbits 

2 

with  such  high  drag  (M^  >  4  deg/day  ),  but  it  is  also  worth  noting  that  the 
orbit  passed  through  15th-order  resonance  with  the  Earth's  gravitational 

field  near  MJD  40600.  Such  a  passage  through  resonance  produces  a  change  in 

•  14  o  • 

inclination  which  for  the  Ariel  3  satellite  was  0.018  :  the  corresponding 

change  for  1969-20B  cannot  be  calculated  because  it  depends  primarily  on  the 

coefficients  and  ^  in  the  Earth's  gravitational  potential. 

*  •  *  0 

However,  the  change  aeems  unlikely  to  exceed  0.003  because  1969- 20B  passed 
through  ’•esonance  so  rapidly.  Nevertheless  an  increase  in  i  of  0.003° 
around  MJD  40600  would  help  to  explain  why  the  last  four  points  in  Fig. 4 
were  above  the  curve. 

5  RESIDUALS  OF  THE  OBSERVATIONS 

The  residuals  of  the  observations,  in  the  form  of  the  observer's  value 
minus  that  given  by  the  fitted  orbit  ct  the  time  specified  by  the  observer, 
have  been  printed  out  using  the  ORES  compute.!  program^  and  have  been  sent 
to  the  individual  observers.  Table  3  gives  the  number  cf  observations  from 
stations  with  5  or  more  accepted  observations,  with  the  rms  residuals  in 
minutes  of  arc,  or  km  for  the  range  measurements. 

Table  3  calls  for  several  comments.  First,  it  should  be  noted  that  the 
directional  residuals  for  the  US  Navy  station  29  are  geocentric  and  should,  be 
multiplied  by  a  factor  of  about  5  to  make  them  comparable  with  the  other 
figures  for  directional  accuracy.  Second,  it  is  possible  that  the  most 
accurate  observations,  namely  those  from  the  Malvern  and  Meudon  cameras,  are 
more  accurate  than  is  suggested  by  the  residuals,  because  this  was  a  high- 
drag  orbit  for  which  a  fit  accurate  to  better  than  2  seconds  of  arc.  though 
often  possible,  cannot  be  guaranteed  if  the  drag  is  irregular.  Third,  the 
proportion  of  rejections  is  rather  high,  and  some  of  these  are  not  necessarily 
bad  observations:  variations  in  drag  on  some  days  may  have  been  so  large  that 
the  polynomial  ;>r  M  was  inadequate  and  good  observations  may  have  been 
rejected,  particularly  if  they  were  isolated  geographically.  Fourth,  the  rms 
residuals  do  not  give  a  realistic  idea  of  the  visual  observer's  capabilities 
because  a  few  observations  made  in  poor  conditions  (semi-daylight,  cloud, 
haze,  etc.)  can  outweigh  a  whole  group  of  good  observations.  For  example, 
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Table  3 


Analysis  of  residuals  for  selected  stations 


Station 

1 

US  Navy 

A 

i. 

US  Navy 

3 

US  Navy 

4 

US  Navy 

5 

US  Navy 

6 

US  Navy 

29 

US  Navy 

414 

Capetown 

2265 

Farnham 

2295 

Leyton 

2303 

Malvern  Hewitt 

Camera 

2304 

Malvern  Radar 

2344 

Thames  Ditton 

2366 

Earl vb urn 

2385 

Stevenage  1 

2403 

Stevenage  3 

2409 

Ribbleton 

2413 

Norwich 

2414 

Bournemouth 

2420 

Willowbrae 

2512 

Mountcastie 

2525 

Crawley  Ridge 

2539 

Dymchurch 

2547 

Luqa 

2573 

Genoa  1 

2577 

Cape  Kinetheodolite 

2586 

Meudon  2 

3101 

Meudon 

3103 

Bordeaux 

8030 

Meudon  Camera 

Others 

TOTAL 

Mo.  of 
observations 


Total  Rejected 


Range 

ktn 


rms  residuals 


Minutes  of  arc 


Dec.  Total 


2.3 

4.7 

2.0 

2.5 

4.3 

8.7 
7.2 

3.1 

2.7 

2.8 
2.0 

3.6 
15.0 

2.5 

5.1 

4.8 

2.9 

0.12 


*Geocentr  Lc 


12 


20 


10  observations  each  with  a  residual  of  1  minute  of  arc  together  with  2 
observations  each  with  a  residual  of  7  minutes  of  arc  give  an  rms  of  3  minutes, 
which  is  a  poor  approximation  to  the  observer's  obvious  capability  of  achiev¬ 
ing  1  minute.  A  better  measure  of  the  observer's  capability  is  the  abc  residual 
(average  in  best  conditions)  the  dfithmetic  mean  of  the  residuals  on  the  best 

70%  of  the  obr  ’jtions.  The  abc  residual  is  usually  about  hal*  the  rms 

.  ,  ,10 

residual  . 

6  CONCLUSIONS 

The  28  orbits  of  1969-2GB  determined  using  the  PR0P5  program  are  of  very 

good  accuracy,  particularly  when  observations  were  available  from  the  Malvern 

Hewitt  camera  or  the  Meudon  camera.  One  of  the  orbits,  No. 14,  with  sd  in 

eccentricity  of  0.00001  and  in  inclination  of  0.0001°,  has  the  lowest  sd  yet 

2 

obtained  for  such  a  high-drag  orbit  (h  =  2.6  deg/day  );  while  orbit  No. 7 
has  sd  in  inclination  of  0.00006°. 

The  values  of  inclination,  after  correction  for  luni-solar  perturbations 
and  odd  zonal  and  tessera1,  harmonics  in  the  geopotential,  have  been  U3ed  to 
determine  zonal  wind  speeds  in  the  upper  atmosphere  at  heights  of  210-250  km, 
averaged  over  latitudes  up  „o  about  25°.  To  fit  theoretical  curves  to  the 
obseived  values  of  inclination,  the  data  has  to  be  split  into  3  parts,  each 
with  a  different  value  for  the  atmospheric  rotational  speed.  The  results  show 
a  clear  distinction  between  che  zonal  wind  at  night  (21  h  to  03  h) ,  which  is 
west-to-east  with  average  speed  140  ±  50  m/s,  and  the  wind  by  day  (04  h  to 
19  h) ,  which  is  east-to-west  with  average  speed  110  ±  50  m/s. 


10  observations  each  with  a  residual  of  1  minute  of  arc  together  with  2 
observations  each  with  a  residual  of  7  minutes  of  arc  give  an  rms  of  3  minutes, 
which  is  a  poor  approximation  to  the  observer's  obvious  capability  of  achiev¬ 
ing  1  minute.  A  better  measure  of  the  observer's  capability  is  the  abc  residual 
(average  in  best  conditions)  the  arithmetic  mean  of  the  residuals  on  the  best 
70%  of  the  observations.  The  abc  residual  is  usually  about  half  the  rms 
residual^. 

6  CONCLUSIONS 

The  28  orbits  of  1969-20B  determined  using  the  PR0P5  program  are  of  very 

good  accuracy,  particularly  when  observations  were  available  from  the  Malvern 

Hewitt  camera  or  the  Meudon  camera.  One  of  the  orbits.  No. 14,  with  sd  in 

eccentricity  of  0.00001  and  in  inclination  of  0,0001°,  has  the  lowest  sd  yet 

2 

obtained  for  such  a  high-drag  orbit  (n  ■  2.6  de?/day  );  while  orbit  No. 7 
has  sd  in  inclination  of  0.00006°. 

The  values  of  inclination,  after  correction  for  luni-solar  perturbations 
and  odd  zonal  and  tesseral  harmonics  in  the  geopotential,  have  been  used  to 
determine  zonal  wind  speeds  in  the  upper  atmosphere  at  heights  of  210-250  km, 
averaged  over  latitudes  up  to  about  25°.  To  fit  theoretical  curves  to  the 
observed  values  of  inclination,  the  data  has  "■<  be  split  into  3  parts,  each 
with  a  different  value  {'•^  the  atmospheric  rotational  speed.  The  results  show 
a  clear  distinction  between  the  zonal  wind  at  night  (21  h  to  03  h) .  which  is 
west-to-east  with  average  speed  140  ±  50  m/s,  and  the  wind  by  day  (04  h  to 
19  h) ,  which  is  east-to-west  with  average  speed  110  ±  50  m/s. 


Table  1  -  ORBITAL  PARAMETERS  FOR  COSMOS  2 ( 


Run 

MJD 

Date 

a 

e 

i 

— 

o 

. 

1 

40316.0 

1969 

7496.878 

0.11919 

48.400C 

199.609 

t 

4 

2 

40328.0 

April  5 
1969 

2 

7468.866 

1 

0.11570 

19 

48.3998 

6 

152.943 

; 

*  3 

40335.0 

April  17 
1969 

2 

7455.106 

4 

0.11411 

22 

&8.395J 

6 

125.488 

; 

*  4 

40341.0 

April  24 
1969 

6 

7444.314 

2 

0.11303 

6 

48.3944 

5 

101.828 

*  5 

40347.0 

April  30 
1969 

2 

7434.111 

8 

0.11213 

7 

48.4014 

6 

78,052 

6 

40354.0 

May  6 

1969 

4 

7421.192 

10 

0.11086 

37 

48.3976 

12 

50.198 

**  7 

40405.0 

May  13 
1969 

3 

7321.433 

6 

0.09902 

25 

48.3848 

7 

201.979 

**  8 

40411.0 

July  3 
1969 

<  1 
7311.948 

5 

0.09782 

1 

48.3835 

1 

176,972 

9 

40420.0 

July  9 
1969 

1 

7299.095 

12 

0.09601 

4 

48.3872 

4 

139.203 

10 

40434.0 

July  18 
1969 

2 

7287.498 

22 

0.09426 

38 

48.3878 

11 

80.214 

11 

40442.0 

Aug  1 

1969 

3 

7279.018 

2 

0.09339 

23 

48.3850 

4 

46.355 

12 

40456 .0 

Aug  9 

1969 

4 

7260.514 

3 

0.09170 

26 

48.3824 

3 

346.754 

13 

40484.0 

Aug  23 
1969 

2 

7210.246 

1 

0.08618 

16 

48.3732 

3 

225.671 

**14 

40492.0 

Sept  20 
1969 

3 

7193.531 

3 

0.08385 

13 

48.3740 

5 

190.574 

15 

40497.0 

Sept  28 
1969 

2 

7180.211 

1 

0.08224 

1 

48.3768 

4 

168.469 

Oct  3 

4 

4 

14 

3 

16 

40503.0 

1969 

7166.513 

0.07995 

48.3719 

141.800 

Oct  9 

2 

4 

20 

4 

17 

40522.0 

1969 

7126.262 

0.07445 

48.3738 

56.312 

Oct  28 

2 

2 

19 

4 

*  18 

40532.0 

1969 

7102.478 

0.07172 

48.3727 

10.687 

Nov  7 

2 

3 

5 

2 

19 

40539 .0 

1969 

7083.106 

0.06966 

48.3684 

338.432 

Nov  14 

3 

2 

15 

2 

20 

40554.0 

1969 

7033.314 

0.06410 

48.3603 

268.202 

Nov  29 

3 

2 

19 

3 

21 

40561.0 

1969 

7009.216 

0.06129 

48.3605 

234.832 

Dec  6 

3 

3 

16 

7 

22 

40574.0 

1969 

6963.248 

0.05524 

48.3529 

171.854 

Dec  19 

2 

1 

13 

3 

Table  continued  on  page  14.  For  key  to  symbols,  see  page  14. 
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FUR  COSMOS  2 >8  ROCKET  C1969-2QB)  WITH  STANDARD  DEVIATION 


c  1 


k  J 

I] 


199.609 

6 

152.943 

6 

125.488 

5 

101.828 


6 

78.052 

12 


201.979 

1 

176.972 

4 

139.203 


346.754 

3 

225.671 

5 

190.574 

4 

168.469 

3 

141.800 


338.432 

2 

268.202 

3 

234.832 

7 

171.854 

3 


0) 

M0 

205.13  , 

22.67 

1 

1 

247.54 

1.10 

2 

3 

272.67 

99.72 

3 

3 

294.10 

100.72 

7 

8 

315.99 

163.56 

7 

9 

340.95 

135.50 

2 

2 

168.64 

163.05 

1 

<  1 

191. z2 

242.74 

1 

1 

'*5.58 

289.87 

3 

4 

279.48 

334.96 

2 

2 

310.36 

234.51 

2 

2 

4.56 

256.18 

1 

1 

113.75 

229 .66 

3 

2 

145.39 

92.78 

1 

1 

165.36 

179.48 

1 

1 

189.64 

81.69 

2 

2 

267.99 

88.08 

2 

3 

324.40 

2 

2 

224.26 

1 

1 

43.22 

56.65 

1 

2 

73.48 

159.57 

3 

3 

4814.414 

2 

4841.521 

2 

4854.929 

6 

4865.490 

2 

4875.508 

4 

4888.244 

3 

4988.478 
<  1 
4998.186 
1 

5011.393 

2 

5023.359 

3 

5032.139 

4 

5051.385 

2 

5104.295 

3 

5122.095 

2 

5136.352 

5 

5151.084 

r 

4. 

5194.782 

2 

5220.893 

3 

5242.323 

3 

5298.083 

3 

5325.426 

3 

5378.239 


2 


1.0736 

11 

0.8690 

15 

0.8916. 

29 

0.8306 


19 

1.0513 

30 

0.7726 

1 

0.8400 

6 

0.5211 


-0.0066 
1 

-0.0089  j  -0.00100 
2  6 


-0.0124 


0.7677 

5 

0.9676 

12 

1.2768 

21 


1.2980 


1.2694 

11 

1.3519 

19 

1.6193 

14 

2.0350 

22 


28 

2.3276 

10 


0.0280 

7 


-0.0048 

1 

-0.0167 


6 


0.0073 


0.0137 

1 

0.0556 

8 

-0.0289 

23 

0.0100 

5 

0.0043 

2 

0.0119 

5 


0.00090 

9 


.91  55  j  8.9 

.01  42 

.84  32  4.0 

.0  1  30  7.0 

.93  32 

.52  19  2.1 

.88  54 

.69  66 

.78  45 


■ 


6 

0.0165 

2 


0.00267 

54 


-0.00067 

6 

0.00553 

24 


0.00178 

37 

-0.00718 

41 

0.00280 
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Fig. I  Values  of  perigee  distance  a  (I- e)  relative  to  Earth  of  radius  R  (=6378*2  km)  *5 


Fig. 3  Values  of  orbital  inclination  L  on  the  28  orbits,  with  sd 


Fig.4  Values  of  inclination  corrected  for  luni-solar,  zonal  and  tesseral 
harmonic  perturbations,  with  fitted  theoretical  curves 


